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evaluation of orbits with incomplete KNO^WiEDGE OF THE 

MATHEMATICAL EXPECTANCY AND THE MATRIX OF COVARIATION OF ERRORS 
D. Ts. Bakhshiyan, R. R. Nazirov, and P. E. El'yasberg 


Examined herein is the problem of selection of the optimal 
algorithm of filtration and the optimal composition of measure- 
ments, assuming that the precise values of the mathematical ex- 
pectancy and the matrix of covariation of errors are unknown. 

In this case, optimization is carried out from the condition of 
attainment of a maximum guaranteed reliability H . of deter- 

g !U’ 

minat.ion of the scalar parameter. H,, ^ . is understood to mean 

ga 1 

the minimum value of reliability, with a given maximum error, in 
a set of possible laws of distribution of the sumn\ary errors 
of the utilized mathematical model of motion and measurements. 
This set is determined by means of the application of some limi- 
tations to the mathematical expectancy and the value of the ele- 
ments of the covariation matrix of errors. The expression for 
, may be utilized for obtaining the guaranteed evaluation of 
the accuracy with, determination of the orbits by a random linear 
rzable algorithm of filtration. It is shown that the problem of 

optimization of H amounts to the solution of some problem of 

gar 

quadratic programming. The optimal algorithm of filtration, ob- 
tained in this case, may be utilized for making some parameters 
more precise {for example, the parameters of the gravitational 
fields) after preliminary determination of the elements of the 
orbit by a simpler method of processing (for example, the method 
of least squares) . 

1. Evaluation of accuracy of Determination of Orbits 

The determination of the orbits amounts to the obtaining of 
an evaluation of q of some vector q={gT7T„,^^ which makes it 
possible to calculate the motion in accordance with the given 
mathematical model [1] . This evaluation is described as the 
function 
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ol; the vector of the measuroments which way be repre- 

sented in the form 


d =d(<i) + ^ , 


( 2 ) 


Where f, is the vector of errors of the initial data, equal to the 
sum of the errors of the measurements and the model (2) [1] . The 
method of least squares (ra.n.k.) is the most widespread algorithm 
of evaluation, with which 


»** 

^ =argmUi [c( - d ( <j !)] '’a 


( 3 ) 


Here, K is the non-negatively determined matrix. The following 
assumptions on the mathematical expectancy of error C and its 
covariation matrix are usually utilized during the derivation of 
the statistical properties of the evaluation of (3) : 


/?(■£, ) = o, z>(4)=ct“a' 




With the assumptions in (4), the evaluation of (3), in many 
cases, is justifiable [1] , i , e. , with an increase in the number of 
utilized measurements, it coincides, in probability, with the true 
value of q. .. of the parameters q:g-^q. with n— +■«>. In reality, 
the assumptions in (4) are never precisely fulfilled, and, with 
sufficiently large n, the evaluation of the accuracy of deter- 
mination of the parameter q, obtained on their basis, proves to 
be unjustifiably optimistic. Depicted in figure 1 is the dependence 
of the root-mean-square value of 0 (E) of the evaluation E of some 
scalar parameter £ on the number n of utilized measurements. In 
this case, the theoretical dependence, which corresponds to the 
assumptions in (4) , is depicted by the dotted line, and the depen- 



dence obtained in practice is depicteid by the solid line. Such 
a situation is inadmissible during the solution of problems of 
navigation or determination of important physical constants, as 
well as during the selection of the optimal strategy of determin- 
ation of an orbit and its correction. 

2, Guaranteed Approach in the Problem of Evaluation of Accuracy 

In connection with what has been set forth, the problem 
arises of obtaining reliable (although less precise) evaluations 
of the accuracy of determination of the orbit parameters. For 
this purpose, the problem of obtaining and optimizing the guaran- 
teed evaluations of accuracy, with incomplete information on the 
statistical characteristics of the errors in the inital data, 
is examined in a number of Studies [1-6] . The present studi- is 
a generalization and further development of the methods of ob- 
taining and utilizing guaranteed evaluations of the accuracy. 

In this case, the results, obtained earlier, may be viewed as 
its partial cases. 

Let the set^Tof possible functions of determination P(5) of 
the error 6 be given 

AY4)e#'. (5) 

and H (or a) is some scalar characteristic, to an increase (or 
decrease) in which corresponds an increase in the accuracy of deter- 
mination of the orbit. We will determine the guaranteed value of 
these characteristics from the inequalities 

( 6 ) 

^ear « i'll " I «gar ^ 

If equalities are achieved in (6) , then we will state that 

and are strict guaranteed characteristics. The charac- 

teristics in (6) are determined for the given algorithm of filtra- 
tion. It is natural to pose the problem of finding of the algorithm 


which is optimal in the sense of an increase (or decrease) in 
the guaranteed (desirably strict) characteristics, i.e. , to 
seek it in a set S of given algorithms of filtration 


"gar 
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In order to calculate the extremums in (6) and (7) , it is 
necessary to prescribe the sets,^and S and the characteristics 
of accuracy. 


3. The Set«.^ 

We will assume that the conditions which determine the set 
^(^contain the following limitations on the mathematical expec- 
tancies E(g^) and the elements of the covariation matrix of 
D(fJ (i,j=r7lT) 


'K(4,)l S M,. , 
o'fj -Vi,- S o,. s D'jj + Vij , i *i ; 0 « ; =5 i>. • + 


( 8 ) 

( 9 ) 


Where M, , , v. .^0 are the elements of the given matrices M, D, 

X X J X J 

V of dimensions Ixn, nxn, nxn. In addition, we will indicate 

the results which take place for che case when, instead of (9) , 

similar limitations are given on the coefficients of correlation 

k. .=D. ./‘^dTTdTT and dispersion: /6 

1 3 3 j ' 1 1 j 3 — 


A* - 
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k ^ < k 
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1/ > 




where k’f'., w. . are the elements of the given matrices K* , W, in 
ij ij 

which k’(‘.=l, w..=0, and i=l,n. 

In addition to conditions (3) -(10), which we will call the 
basic conditions, the set ^^may be determined by additional limi- 
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tations aa well, which can be placed on the function of distri- 
bution of the errors (for example, the form of this function may 
be prescribed) , 

We will give examples of the representation of the set^^ 
using (8) -(10). With M*0, V«*0, we obtain the classic assumptions 
in (4), on the basis of which we construct the evaluations of 
the method of least squares [7]. If 14^0, D*»V»0, then the set 
^is characterized by the condit:‘ons (the "skirting" 

scheme) [2]. The cases ki«kj^<?k 2 , and lk^j|*k<fl, 

are examined in [1,5,6]. 


hexr Guarant 


In the capacity of the characteristic of accuracy H, we will 
take the probability of the fact that the error 




(n) 


of determination of some scalar parameter il(q) does not exceed 
the given magnitude of a>0, according to the modulus, i.e. , 


H = P ( \ SC \ ^a, ) . 


(I?) 


If H is given with reliability, then one can view a=a(H) as 
the characteristic of accuracy. The solution of problems (6) , 

(7) for the introduced characteristics H and a, with the conditions /l^ 
given below, are equivalent [1] . Therefore, we will subsequently 
make use of the characteristic H . 

We will calculate H, , . for the algorithms of filtration (1) , 

gai' 

which satisfy the following conditions. 

1. If 5=0, then 6«.=0 for all q from the region of its possible 
values. According to (1) , (2) , and (11) , this condition may be 
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written in the form 




If (13) is fulfilled, then we will state that the algorithm of fil- 
tration (1) possesses the property of unbiasedness. 


2. The function £-tq(d)] is linearizable at some prior known 
point , i . e , , 


= <■(<,„) +.V[rf 


(T4) 


Where ) is a matrix of dimensions 1 x n. Hence, taking 

(13) into account, it is not difficult to obtain 




(T5) 


We will find H for two cases, 
gar 

A. In addition to conditions (8) , (9) , there is no additional 
information on F(F,), which makes it possible to find the law of 
distribution of the error fia. 


B. The distribution Si may be considered normal. (Such an assump- 
tion corresponds to the truth if, for example, the errors 5 satisfy 
the conditions of the central limit theorem) . 


For case A, utilizing Chebyshev's inequality, we obtain the 
following from (6) for the strict guaranteed characteristic: 


H 


•gar 


where 




= '"gar *^gar' 


/8 
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are strict guaranteed characteristics. 


In case B, the reliability of (12) is equal to [1] 


where 


Jl = If (a., rr,,l>) = <P ~ ^ > 

7Tl ^\K{t>f)\, D^p(t), 


(To) 


We will assume that <%>m, since, in the opposite case, 

H ((x,m,D) <0 ,5 , which is unacceptable in practical problems of 
the determination of orbits. Then, the derivatives of function 
(18) , according to m and D, are negative, in virtue of which wo 
obtain 


Ml « n H (a, 


m 


i> ) “ R, 




■^sur) 




v;here m„„„ , D„„ , are found from (17). 
gai gai 

Thus, in both cases A and B, for the calculation of H , 
it is necessary to find the strict guaranteed characteristics (17) . 
According to (8), 


^gar “ -f-/M f 


( 20 ) 


where X = ( ]xx [ , | x I ) . During the calculation of D , it should 

n 8^^ 

be taken into account that the acceptable matrices D(C) should 
satisfy not only conditions (9) , but should also be non-negatively 
determined. In this connection 


where the right-hand portion in (21) is the maximum of D(£), with 
the conditions in (9) . 



For D and V, the necessary and sufficient conditions may be 
obtained so that in (21 ) , with any X, the equality would be ful- ^ 

filled. It follows from them, specif icaHy, that the conditions 
of non-negative definiteness of the matrices D*+Z and V-Z are 
sufficient with some value of the diagonal matrix Z, vjhile the 
non-negative definiteness of the matrix D*+V is necessary. 

Now, let the limitations on the covariations be given in the 
form of (10), and, under these conditions, either w. .>jkt.l, i?*j, 

A J X J 

or In this case, we will adopt the designations A* 

diag [Ai , , AD*»*6K* . Then, one can show that the maximum of 

D(l), with the conditions in (10), may also be written in the form 
of the right-hand portion of (21) , This makes it possible to 
subsequently view only that case when tae conditions in (9) are 
given, since similar results will occur with limitations in the 
form of (10) . 

Sufficient conditions for the implementation of the equality 
in (21) are fulfilled for all of the examples enumerated in 
paragraph 3. Subsequently, we will also assume that these con- 
ditions are fulfilled, and there is an equal sign in (21) . 

5. Optimization of the Guaranteed Reliability 

Relationships (16) -(21) give strict guaranteed values of the 
characteristic H for cases A and B from paragraph 4. Therefore, 
in both variants, problem (7) is reduced to maxi-mini optimization 

6 

tnAX II = Tnajc rnin H. (2?) 

X X 

Here, the maximum is taken according to the line X, which, 
according to (14) , determines the algorithm of filtration and 
satisfies the condition of unbiasedness (13) of this algorithm. 

We will assume that the functions d(q) and <t(q), similar to 
(14) , are linear /IQ 
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where h an^ C are known matrices of dimeni<i.ons n x m and I x n. 
Prom (11), (14), and (23), we obtain that the condition of un- 
biasedness (13) takes on the form 


XA =c . 


( 2 “) 


For case A, problem (22), according to (3 <j) and the obtained 
strict characteristics, amounts to finding 


U-C 


gar 


^in (X/)">^+.V+(/MM''4| )X+} . 



Similar to that, as was done in [5], minimization in (25), by 

ir«PCiii23 Qi4 sui.jaua.uUu4.On Qj vU4« 4.ciAjJ.\=a 


* x; - 4c; , I a:,. I «.t; 4 a:/ , = f 


a 


( 26 ) 


is brought to the problem of quadratic programming, the solution 
of which satisfies the conditions x|x?»0, which are necessary for 
the correctness of the substitution (26) . 


We will now switch to case B. In order to find the maximum 
in (22) , we will examine first the auxiliary problem 


max M 


(27) 


with the conditions in (24) and the additional limitation 


m 


gar 




(26) 


where m is determined from (20) , and ).i=const, 
gar 


If the solution of this problem will be known for all of the 



accQpfcabl© li 


/II 


H-f l"»o , W/,1 


>r»o* >n rwi -cot * 

)(A^c 



then the maximum (%Z\ is found by means of unidimensional optimi- 
zation according to u in the interval in (29) . Since the deriva- 
tive of according to D is negative/ the maximization in 
(27) amounts to the finding of 


^(H)= v«in po V+;c* \/Xl\ 


with the conditions (24) and (28), Thus, problem (22), in case B, 
is reduced to unidimensional maximization, according to w, of 
the function H(n,u,f(w)) in the interval (29) (see fig. 2), The 
finding of the function i{u) in this segment, after the substi- 
tution (26) , amounts to the problem of quadratic programming with 
the parameter p in the right-hand portion of condition (28) . The 
last problem is of the very same order of crmplexity as the prob- 
lem of quadratic programming (8,9). During its solution, the 
interval (29) is subdivided into a finite number of segments, 
inside of each of which the solution is acbreved on one and the 
same basis, and f(p) is an arc of a parabola. One can show 
that, in each of these segments, there is no more than one critical 
point, of the function Il(a,p,f (p) ) , and this point is found by 
means of a search of the root of some increasing function f (p) 

(see fig,. 2) . A simple graphic method is developed along with 
the described analytical method for finding the optimal p. 

Thus, for the examined practically important cases A and B, 
the solution of the problem of optimization of the guaranteed 
reliability may be obtained using the well-developed algorithms 
of quadratic programming. 

6. Determination of the Radial Velocity According to Range /12 

Measurements 

In the capacity of an example, we will examine the case when 
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m roust dotemine bh® radial velocity v of a cosmic object at 
a moment in time b»0, according to measurements of the ranges 
from a measuring point at the moment t^. In this case, it 
is assumed 'that some prior dependence a(t) of the radial velocity 
on time is known. Then, relationship (2) has the form 


^ t t * 

d; » r + S «(T) rft + 1 i , 

0 0 


( 51 ) 


whore r is the radius with t»0, and is the sum of the error 
of the measurements and the error evoked by the inaccuracy of 
knowledge of the radial acceleration, which we will assume is 
a non-random magnitude, not exceeding the given number w. We 
will also assume that the mathematical expectancy of errors of 
the measurements and the coefficients of correlation between 
them, for all errors, are limited according to the modulus 
by the given magnitudes M>0 and K^O, and the dispersions do not 
exceed <r^. Then, for the summary error we obtain 


We will assume that, at each moment t^, repetition of measurements 
is allowed, so that their total number would not exceed some n. 


One can show that the optimal solution of problem (22) , for 

cases K and B, consists of the conduct, according to n/2, of 

measurements at some moments ti«T and tg“-T. in this case, 

max may be represented as a function of the dimensionless 

parameters a-a/0, M=M/y, where B^s=minE((SS,^) „ ^ and t (l-k)/n+k]/o*'’ . 

X gar 

Given in figure 3 are the graphs of the dependences max H (a) , 
found for values of M«0,l,2; in this case, the solid lines corres- 
pond to case A, and the dotted lines to case B. It is evident from 
figure 3 that the sv/itch from case A to case B involves a slight 
decrease in the maximum error a (no greater than 17% with equal 
Hj,^,s0.999). 


Calculations were carried out for the examined example, which 
showed that, with filtration of the method of least squares, with 
equal weights in some segment obtained 

only then differs little from its optimal value, when the val»>e of 
is correctly selected. With other the guaranteed 

reliability decreases sharply with evaluation of the method of 
least squares. 

7. Practical Utilization of the Optimal Algorithms of Filtration 

The method, proposed in the present study, for constructing 
the optimal algorithm of filtration makes it possible to increase 
the guaranteed accuracy of determination of the orbit, as compared 
with the results obtained using the classical algorithms, based 
on the given values of the mathematical expectancy and the co- 
variation matrix of errors. However, it possesses a number of 
shortcomings. The basic of these are; 

-considerably greater (as compared with classical algorithms) 
labor consumption of calculations; 

-dependence of algorithm on selection of parameter being 
evaluated ; 

-optimalness of algorithm only with linear posing of the 
problem, leading to the necessity of solution of nonlinear prob- /14 
lems by the method of iterations. 

The indicated circumstances evidently limit the utilization 
of the proposed optimal algorithm of filtration for the solution 
of operational problems of determination of the param<=^ters of 
motion of a space vehicle. However, with secondary processing 
of the data of trajectory measuremients , in order to determine the 
parameters which have independent scientific value (elements of 
orbits of natural heavenly bodies, parameters of gravitational 
fields of these bodies, and so on), the use of the indicated optimal 
algorithm of filtration may prove advisable. In this case, the 
follov^ing order of processing of the measurements may be proposed. 

1. The prelimi.naiy orbit, which is utilized as the reference 
orbit for construction of a linearized mathematical model (2) , is 
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determined by one of the claasical methods. 


2. According to the results of the study of the utilized 
measurement system and mathematical model, and also on the basis 
of the correlation analysis of the discrepancies of the system of 
conditional equations, obtained in the process of preliminary 
determination of the orbit# one can find the possible limits of 
change in the mathematical expectancy and the covariation matrix 
of errors. According to these data, conditions are selected which 
determine the set,0^ and the exclusion of anomalous measurements 
is carried out., 

3. fhe value of each of the interesting parameters Z is 
made more precise, using the corresponding optimal algorithm of 
filtration, obtained from the solution of the problem of quadratic 
programming « 

4. The guaranteed characteristics of accuracy of the obtained 
evaluations of i are determined according to formulas (20) and (21) . 
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FIG. 3 

Dependence of reliability of max H on the dimension- 
less maximum error a with various values of the parameter M 
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